'Афуапсеа Епэтеетия Кезеагсй 2021. Т. 21, № 2. С. 191-199. 155М 2687—1653 





ИНФОРМАТИКА, ВЫЧИСЛИТЕЛЬНАЯ 
ТЕХНИКА И УПРАВЛЕНИЕ 
1УМЕОВМАТТОМ ТЕСНМОГОСУ, СОМРОТЕВ 
ЭСТЕМСЕ, АХО МАМАСЕМЕМТ 











УДК 620.179 Брз://401.0г2/10.23947/2687-1653-2021-2 


Алгоритм управления упруго-вязко-пластичной моделью для 
исследования процессов ударного взаимодействия тел 
В. Л. Лапшин, Е. В. Зеньков"? [=] ИН = 


1 ФГБОУ ВО «Иркутский национальный исследовательский технический университет» (г. Иркутск, Российская Федерация) 
? ФГБОУ ВО «Иркутский государственный университет путей сообщения» (г. Иркутск, Российская Федерация) 





Введение. В инженерной практике важное значение имеют динамические процессы, с помощью которых опи- 
сывают и изучают механику взаимодействия деталей машин и элементов конструкций. Эти динамические про- 
цессы являются причиной возникновения больших деформаций, приводящих к разрушению. Целью исследова- 
ния являлось создание алгоритма моделирования ударного процесса за счет управляемого преобразования ме- 
ханореологической модели ударного процесса из упруго-вязкой в упруго-вязко-пластичную. 

Материалы и методы. Предложены дифференциальные уравнения движения модели. Рассмотрены условия 
преобразования модели при переходе от упругих к пластическим деформациям, от этапа нагружения модели к 
этапу разгрузки при ударном взаимодействии с поверхностью. При расчете деформаций используется допуще- 
ние о том, что упругие и пластические деформации возникают одновременно с самого начала удара. Подробно 
рассмотрена методика функционирования модели, составлен алгоритм работы модели, подробно описана логи- 
ка его функционирования. 

Результаты исследования. Для исследования ударных процессов была разработана механореологическая 
упруго-вязко-пластичная модель. Важным параметром модели является усилие, соответствующее началу обра- 
зования пластических деформаций. В результате исследований был создан более совершенный алгоритм и раз- 
работана новая компьютерная программа для исследования ударного процесса с помощью упруго-вязко- 
пластичной модели с регулируемым упруго-пластическим преобразованием. 

Обсуждение и заключения. Полученные результаты могут быть использованы для повышения точности и до- 
стоверности моделирования ударных процессов с целью дальнейшего развития методик и способов определе- 
ния физико-механических характеристик материалов ударными методами. Знание механических характеристик 
материалов необходимо при решении различных исследовательских задач путем математического моделирова- 
ния вибрационных и ударных процессов. При этом важной задачей является адаптация расчетной модели к ре- 
альному ударному процессу, для чего необходима разработка соответствующих способов и методик. 
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ние ударного процесса. 
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Введение. В инженерной практике важное значение имеют динамические процессы (ударные, 
вибрационные). С их помощью описывают и изучают механику взаимодействия деталей машин и элементов 
конструкций. При этом большое внимание уделяют ударным процессам, так как при этом тела подвергаются 
действию больших динамических нагрузок. Они могут являться причиной больших деформаций и приводить к 
разрушению оборудования. Данные обстоятельства объясняют повышенное внимание к моделированию и 
исследованию динамических процессов со стороны ученых и инженеров [1-14]. Особый практический интерес 
имеет развитие экспресс-методик и способов определения физико-механических характеристик материалов 
ударными методами. При этом важной задачей является адаптация расчетной модели к реальному ударному 
процессу. В связи с этим целью исследования является создание более точного алгоритма моделирования 
ударного процесса за счет управляемого преобразования механореологической модели ударного процесса из 
упруго-вязкой в упруго-вязко-пластичную. Разработанный алгоритм позволяет на этапе нагружения 
осуществлять преобразование модели на основе заданных условий, что обеспечивает более качественное и 
достоверное моделирование ударного процесса с учетом реальной пластичности материала. Предыдущая 
разработка не позволяла осуществить данное преобразование и модель при ударном взаимодействии являлась 
упруго-вязко-пластичной с самого начала нагружения. 

Материалы и методы. Одним из направлений научных исследований в Иркутском национальном 
исследовательском техническом университете является моделирование и изучение процессов динамического 
взаимодействия твердых тел. Для изучения динамики протекания данных процессов с целью решения 
практических задач разработана и исследована упруго-вязко-пластичная модель [15]. Модель построена с 
использованием реологических и инерционных элементов т! и то, объединенных в два блока (рис. 1). Один из 
блоков является упруго-вязким К!-С. Он предназначен для моделирования упругих деформаций тела. Они 
возникают под действием внешней нагрузки и исчезают при ее отсутствии. Упруго-пластический блок К>—№ 
моделирует пластические деформации, остающиеся после снятия нагрузки. Масса тела учитывается 
инерционным элементом ти. Масса т› введена для того, чтобы описывать динамику движения упруго-вязко- 
пластичной модели с помощью двух дифференциальных уравнений второго порядка. Чтобы исключить 
влияние 712 на динамику процесса, ее величина принималась малой (т < 0,001 т/). Как показали расчеты, при 
таком условии можно принять №1= (4). 

Рассмотрим ударный процесс на примере сферического тела. При нагружении тела в нем возникают 
внутренние силы, препятствующие изменению формы и размера тела, то есть препятствующие (2) его 
деформации. Одной из них является сила сопротивления упругой деформации. При этом могут возникать 
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дополнительные сопротивления деформации, которые называются вязкими. Они связаны со скоростью 
деформации тела. На модели данные сопротивления деформации тела описываются выражениями: 
№ р Е, р, 12 


Е =Со-у,"-,", (1) 


В =К(у,- У, И 


где у, у, у, у, — динамические характеристики модели, описывающие перемещение и скорость инерционных 


элементов т и 12; К! — параметр модели, характеризующий упругие свойства материала (упруго-вязкого 
блока модели); С — параметр модели, характеризующий вязкие (диссипативные) свойства материала в том же 
блоке модели (3). 
Силы сопротивления деформации, возникающие в упруго-пластическом блоке модели, описываются 
следующими выражениями: 
М, = Е, +Е,,,М, 2 №; 


п2 
Е (2) 
из 

Ра = У У, + бы 
Здесь Е —Щ_ усилие нагружения модели, при котором начинает работать упруго-пластический блок модели, 
учитывающий пластические деформации материала; К› — параметр модели, характеризующий упругие 
свойства материала при пластической деформации (упруго-пластического блока модели); № — параметр 
модели, характеризующий пластичность материала в том же блоке модели (3). 

















Рис. 1. Упруго-вязко-пластичная механореологическая модель 


Переменные у,у описывают величину и скорость полной (суммарной) деформации. Переменные 


уу, описывают величину и скорость пластической деформации. Следовательно, скорость и величина 


упругой деформации определяется как у-—у,; у/-у,. 

Реологические элементы модели (упругий, вязкий, пластичный) имеют нелинейные характеристики, то 
есть упругие, вязкие и пластические сопротивления деформации пропорциональны скорости и величине 
деформации в заданной степени (51, 52, И1, 2, из). В качестве рекомендации можно отметить, что для упругих 
сопротивлений при ударе сферического тела следует использовать и1=и›=3/2 [2—4], для вязких и пластических 
сопротивлений можно принять значение равное единице [2, 3]. 

Результаты исследования. Рассмотрим работу модели. На начальном этапе использования модели в 
научных исследованиях использовался упрощенный вариант ее функционирования. Принималось условие, что 
при деформации тела упругие и пластические деформации возникают одновременно с самого начала 
нагружения. Это позволило разработать и опробовать на практике более простой алгоритм функционирования 
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модели. На его основе была разработана исследовательская компьютерная программа', выполнен комплекс 
научных исследований ударного процесса на модели”, в том числе представленных в работе [16]. 

Однако данное упрощение ограничивает возможности модели и не позволяет в полном объеме 
выполнять исследования ударных процессов. Поэтому в результате дальнейших исследований был создан 
более совершенный алгоритм и разработана новая компьютерная программа для исследования ударного 
процесса с помощью  упруго-вязко-пластичной модели с регулируемым — упруго-пластическим 
преобразованием". Рассмотрим этот алгоритм функционирования модели более подробно. 

Разработанный алгоритм обеспечивает следующую последовательность функционирования блоков 
модели. На начальном этапе ударного процесса в работе участвует только упруго-вязкий блок модели К!-С. Он 
позволяет рассчитать упругие деформации тела на этапе нагружения (рис. 2). При этом упруго-пластический 
блок модели не функционирует, так как до определенного момента пластические деформации в теле могут 
отсутствовать. Появление пластических деформаций соответствует заданному усилию ЁР5г. Когда расчетное 
усилие нагружения №! достигнет заданного значения (М=ЁР5т), упруго-вязкая модель преобразуется в упруго- 
вязко-пластичную. Начинает функционировать упруго-пластический блок модели. Когда усилие нагружения 
достигает своего максимального значения М№,., рост деформации прекращается и этап нагружения 
заканчивается. За счет упругости материала тело начинает восстанавливать свою первоначальную форму и 
размеры, наступает этап разгрузки модели. В ходе этого процесса исчезают упругие деформации, при этом 
возникшие пластические деформации остаются в теле и, соответственно, на модели. Таким образом на этапе 
разгрузки модели работает только упруго-вязкий блок. 































































































Рис. 2. Схема преобразования модели 


Этап контактного взаимодействия заканчивается, когда сила нормальной реакции №, определяемая 
уравнениями (1) и (2), становится равной нулю. Тело отскакивает от поверхности на высоту й. 
Процесс нагружения модели на различных этапах деформирования описывается с помощью 


дифференциальных уравнений, которые приводятся ниже: 
51 52 











т, У +С, У У, +Ку”" =-—т&, (3) 
т у+С(у- у," У," +К, 0 -У,)" =-тв, (4) 
т, у +Ку" + Ху, + С (уу, +Кб,-у)" =-т,в+Е,. (5) 


С помощью уравнения (3) рассчитываются упругие деформации модели на этапах нагружения и 
разгрузки. С помощью уравнений (4) и (5) рассчитываются упругие и пластические деформации модели при ее 
нагружении на этапе пластического деформирования материала. 

Сила нормальной реакции на этапе упругой деформации определяется по формуле: 


.м 52 


№ =Су, у, +Ку". (6) 


Сила нормальной реакции на этапе пластической деформации определяется по формуле: 
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п3З 


№, = Ку" + у, +Ру. (7) 
Вычисление производится с заданным временным шагом Е. Время удара определяется по формуле: 
Рей (8) 


где 7 — количество шагов. 
Условие перехода от упругой деформации к пластической имеет вид: 


№ > Ех (9) 
Условие окончания этапа нагружения и начала этапа разгрузки модели: 
м, >М, (10) 
Условие окончания этапа контактного взаимодействия модели: 
№ =0 (11) 


Рассмотрим алгоритм исследования ударного процесса (рис. 3). 
1. Блок 1 формирует исходные данные исследуемого процесса. 
2. Блок 2 вводит исходные данные процесса: начальное время удара { = 0, начальную деформацию тела 


(модели) у = 0, начальную скорость деформации у = у, ‚ начальную силу удара М = 0. 

3. Блок 3 выполняет вычисление значений коэффициентов модели Ку, К», №, С. 

4. Блок 4 выполняет вычисление времени удара и фиксирует количество шагов расчета на этапе 
нагружения при упругой деформации. 

5. Блок 5 выполняет вычисление параметров процесса ударного взаимодействия тел с помощью 


упруго-вязкого блока модели. Рассчитывается упругая деформация у! и скорость деформации у, (уравнение 


(3)). Определяется сила ударного взаимодействия М, (уравнение (6)). 

6. Блоки 6 и 16 осуществляют проверку по условию (11). При невыполнении условия расчет 
продолжается на следующем шаге этапа ударного взаимодействия. Если условие выполняется, то расчет 
завершается и выполняется переход к этапу отскока модели. 

7. Блоки 7 и 12 контролируют момент завершения процесса нагружения модели. Выполнение условия 
(10) означает, что этап нагружения модели закончился, сила контактного взаимодействия №; начинает 
уменьшаться. Выполняется расчет на этапе разгрузки модели. Вычисление передается блоку 4 или 14. При 
невыполнении условия (10) сила контактного взаимодействия №; продолжает увеличиваться. Этап нагружения 
модели продолжается и вычисление производят блоки 8 или 13. 

8. Блок 8 фиксирует наибольшую величину силы контактного взаимодействия №„„, соответствующее 
ей время в процессе удара Т»их и передает управление блоку 9. 

9. Блок 9 определяет окончание этапа упругой деформации при нагружении модели. Если условие (9) 
выполняется, значит сила нормальной реакции модели № превышает заданное усилие ЁЕ„, соответствующее 
началу образования пластических деформаций. Этап упругой деформации завершается, начинается этап 
упруго-пластической деформации и управление передается блоку 10. Если условие (9) не выполняется, значит 
сила нормальной реакции модели №; меньше заданного значения Р„, и пластические деформации отсутствуют. 
Вычисление передается блоку 4 для выполнения расчета на следующем шаге этапа упругой деформации. 

10. Блок 10 производит расчет текущего времени и номера шага вычисления на этапе упруго- 
пластической деформации. 

11. Блок 11 выполняет вычисление параметров процесса ударного взаимодействия тел на этапе 


пластического деформирования. Рассчитываются скорости и величины деформаций модели у,у,у,,У, 


(уравнения (4) и (5)), определяется сила контактного взаимодействия №, (уравнение (7)). Далее управление 
передается блоку 12. 

12. Блок 12 определяет окончание этапа нагружения модели. Если условие (10) выполняется, значит 
сила нормальной реакции №, уменьшается, этап нагружения модели закончился и выполняется переход к этапу 
разгрузки модели. Управление передается блоку 14. Если условие (10) не выполняется, значит сила нормальной 
реакции №; продолжает увеличиваться, идет этап нагружения модели и управление передается блоку 13. 

13. Блок 13 фиксирует наибольшую величину силы контактного взаимодействия № и 
соответствующее ей время в процессе удара Т»„„х и передает управление блоку 10 для расчета на следующем 
шаге этапа упруго-пластической деформации. 

14. Блок 14 производит расчет текущего времени и номера шага вычисления на этапе упругой 
деформации при разгрузке модели. Управление передается блоку 15. 
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15. Блок 15 выполняет расчет параметров процесса ударного взаимодействия тел с использованием 


упруго-вязкой модели при разгрузке: упругой деформации у1, скорости деформации у, по уравнению (3), силы 
нормальной реакции М; по уравнению (6). Управление передается блоку 16. 

16. Блок 16 осуществляет проверку условия (11). При невыполнении условия этап ударного 
взаимодействия продолжается и управление передается блоку 14. При выполнении условия расчет на этапе 
ударного взаимодействия завершается и осуществляется переход к этапу отскока модели. Управление 


передается блоку 17. 
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Рис. 3. Схема алгоритма 
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17. Блок 17 определяет высоту отскока тела й, время полета при отскоке Тк и передает управление 
блоку 18. 

18. Блок 18 формирует и выводит основные характеристики ударного процесса: время ударного 
взаимодействия тел (модели) 7 (уравнение (8)); силу ударного взаимодействия тел №»; время в процессе удара 
Тута» соответствующее №„„х; высоту отскока тела (модели) й. 

На базе разработанного алгоритма была создана исследовательская компьютерная программа“. С 
использованием авторской программы были выполнены исследования процесса ударного взаимодействия 
тел [17]. Изучалось влияние различных факторов на динамику ударного процесса. В частности, на модели 
изучалось влияние механических свойств материала (упругости, вязкости, пластичности) на время и силу 
ударного взаимодействия сферического тела с поверхностью, величину и скорость деформации, высоту отскока 
после удара. При этом рассматривались различные варианты упруго-пластического преобразования 
исследовательской модели, отражающие переход из упругого в упруго-пластическое состояния материала при 
нагружении. 

Сравним работу упруго-вязко-пластичной модели по предыдущему (график 1, рис. 4) и новому 
(график 2, рис. 4) алгоритму. 
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Рис. 4. Изменение силы ударного взаимодействия тел № 


Новый алгоритм позволяет при моделировании ударного процесса взаимодействия тел осуществлять 
переход от упруго-вязкой модели к упруго-вязко-пластичной модели по заданному усилию ударного взаимо- 
действия Ест. Этот переход на модели соответствует началу образования вматериале пластических деформаций, 
что позволяет описывать ударный процесс более достоверно и точно. Хорошо видно, что закономерность изме- 
нения силы нормальной реакции, рассчитанная по новому алгоритму (график 2, рис. 4), существенно отличает- 
ся от аналогичной закономерности (график 1, рис. 4), полученной на предыдущем варианте алгоритма. Отлича- 
ется время удара Тк и время от начала удара до максимума силы ударного взаимодействия Туна» (рис. 4). Может 
существенно изменяться и величина №». Таким образом, разработанный алгоритм позволяет более качествен- 
но решать вопросы адаптации упруго-вязко-пластичной модели к реальным процессам ударного взаимодей- 
ствия тел с учетом их механических свойств. 

Заключение. В результате исследований был создан более совершенный алгоритм и разработана новая 
компьютерная программа для исследования ударного процесса с помощью упруго-вязко-пластичной модели с 
регулируемым упруго-пластическим преобразованием. Представлены дифференциальные уравнения движения 
модели. Рассмотрены условия преобразования модели при переходе от упругих к пластическим деформациям, 
от этапа нагружения модели к этапу разгрузки при ударном взаимодействии с поверхностью. Подробно 
рассмотрена методика функционирования модели, составлен алгоритм работы модели, подробно описана 
логика его функционирования. 


* Удар упруго-вязко-пластичной модели с регулируемым упруго-пластическим преобразованием: св-во о рег. прогр. для ЭВМ. 
В. Л. Лапшин, Е. В. Зеньков. № 2019618137; дата регистрации 26.06.2019 г. 
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Преимуществом нового алгоритма является то, что он позволяет более точно описать процесс в момент 
начала образования пластических деформаций в материале при его нагружении. Момент начала пластических 
деформаций на модели может быть точно задан и полностью соответствовать фактическому поведению 
материала. Для этого, например, можно использовать данные экспериментальных исследований ударного 
процесса. В предыдущих версиях алгоритма и программы такая возможность отсутствовала. 

Так, например, для исследования ударного процесса взаимодействия сферического тела с 
поверхностью был разработан алгоритм адаптации упруго-вязко-пластичной модели к условиям ударного 
процесса взаимодействия тел‘. 

В данном алгоритме было реализовано условие, при котором упругие и пластические деформации 
возникают одновременно с самого начала ударного взаимодействия. Однако данное упрощение реального 
процесса ограничивает возможности моделирования, особенно при изучении влияния пластичности материала 
на динамику ударного процесса. 

Усовершенствованный алгоритм позволяет реализовать возможности упруго-вязко-пластичной модели 
в полном объеме. Упруго-пластическое преобразование модели может быть выполнено по заданному плану, 
что обеспечит более качественную адаптацию модели к реальному ударному процессу. При этом фактические 
пластические деформации могут быть смоделированы более достоверно. 

Таким образом, новый параметр модели Ест (усилие нагружения, соответствующее появлению в 
материале пластических деформаций), введенный с помощью разработанного алгоритма, расширяет 
возможности адаптации модели к реальным ударным процессам и повышает точность модели при описании 
пластических деформаций. Работоспособность алгоритма была подтверждена в ходе проведения исследований 
ударного процесса с помощью модели [18]. 

Практическое применение разработанной программы может существенно повысить точность и 
достоверность моделирования ударных процессов. В частности, результаты исследований ударных процессов 
могут являться базой для совершенствования методик и способов оценки механических свойств 
конструкционных материалов, основанных на ударном принципе действия. Так, на основе предыдущей версии 
алгоритма был разработан новый способ определения модуля упругости материала, основанный на ударном 
взаимодействии тел [16]. Знание механических свойств материалов необходимо при выполнении исследований 
для формирования математических моделей процессов вибрационного и ударного взаимодействия 
механических систем. При этом важной задачей является адаптация расчетной модели к реальному ударному 
процессу, для чего необходима разработка соответствующих способов и методик. Решение данных задач 
возможно на основе полученных результатов исследований. 
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